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A complete solution describing a binary system constituted by two unequal counter-rotating 
. black holes with a massless strut inbetween is presented. It is expressed in terms of four arbitrary 

parameters: the half length of the two rods representing the black hole horizons cti and 02, the total 
' mass M, and the relative distance R between the centers of the horizons. The explicit form of this 

solution in terms of physical parameters, i.e., the Komar masses Mi and M2, the Komar angular 
• momenta per unit mass ai and 02, having ai and 02 opposite signs, and the coordinate distance 

' R, led us to a 4-parameter subclass of solutions in which a set of five physical parameters satisfy a 

, simple algebraic relation. Moreover, the interaction force, provided by the strut, between the black 

■ holes results to be of same form as it is for the static double-Schwarzschild case. 

(N 
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I. INTRODUCTION 



The equilibrium configuration of the famous double-Kerr-NUT solution [l| has been extensively studied in the 
last three decades. Applying regularity conditions on the symmetry axis, the balance equations were first derived 
, by Kihara and Tomimatsu [21, Hoenselaers solved them analytically for the case of subextreme sources, i.e., non- 
' degenerate black holes, where its analysis revealed the existence of ring singularities outside the horizons, outside of 
00 ! the symmetry axis due principally to the fact that at least one of the Komar masses is negative [Ij. 

■ Recently Neugebauer and Hennig [6^ have shown the non-existence of regular solutions describing equilibrium 
[ configurations between two black holes, by using the analytical solution presented by Manko et al. 7]. Additionally, 
. if the parameters do not satisfy the regularity conditions on the symmetry axis, there arise two kind of singularities 
CO ' on the axis, known after Bonnor (sj as torsion singularity and stress singularity. The first one generates a region with 
closed timelike curves due to the presence of NUT sources, which lead to finite and semi-infinite singularities along 
the axis, breaking the asymptotic flatness of the solution [9|, idj. 

The second one represents a strut, a conical singularity ll|, which helps us to understand the interaction force 
between the two bodies by means of the gravitational attraction and a spin-spin interaction. 



In the aforementioned equilibrium problem, in the absence of a strut, one always starts by solving the balance 
condition, then the corresponding algebraic variables are substituted into the axis condition. Nevertheless, one might 
choose the opposite way and first solve the equations for avoiding the NUT sources, with the purpose of calculating 
the massless strut and determine the interaction force between the two black holes. This last approach is more general 
and complicated to analyze with respect to the equilibrium situation. Nowadays this 5-parameter subclass of the 
well-known double-Kerr-NUT solution [T] and its analytical representation in terms of physical parameters remain 
still as an open problem. 

One of the first attempts to describe the physical properties of two rotating black holes, was made by Varzugin p^ . 
He solved the corresponding Riemann-Hilbert problem, in which the irreducible mass ai is defined as the half length 
of the rod representing the event horizon of the i-th. black hole located on the symmetry axis. First, the axis condition 
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for the system of two black holes separated by a niassless strut is formulated. The massless strut is introduced in 
order to avoid the falling of the black holes into each other, due to the gravitational attraction. 

Since, by means the Smarr mass formula [l^, the irreducible mass is related to the respective surface gravity 
and area of the black hole horizon, one could express it in terms of physical Komar parameters. Varzugin found 
an analytical solution describing a binary system constituted by identical counter-rotating black holes, where the 
interaction force has the same form as the one for the static Schwarzschild case. The corresponding unique a is 
described by only three parameters, i.e., the Komar's mass m and angular momentum per unit mass a for the upper 
black hole jS]. The lower black hole has parameters m and —a, and R is the coordinate distance between both 
constituents. 

Later on Manko et al. l3] constructed the simplest configuration of two counter-rotating black holes. A 3- 
parametric solution using the explicit form of a. They wrote explicit expressions for the Ernst potential on the 
syinmetry axis and for the complete metric outside the symmetry axis. This solution is equatorially antisymmetric 
[15| . where the axis condition is straightforward fulfilled and the total angular momentum of the system vanishes, i.e., 
J = 0. 

In the framework of the case of two identical bodies, Bonnor 0, [l^ advanced two additional conditions to be 
satisfied in order to remove the additional contribution provided by the massless spinning rods outside the sources: 

0^ + 11 = 0, ^^)a, + a, = ,. (1) 

The first condition of Eqs.([T]) avoids the semi-infinite massless spinning rods, located in the upper and lower parts 
of the symmetry axis, while the second one, avoids the massless spinning rod of finite length between the two bodies. 



On the other hand, it is worthwhile to stress the fact that the most satisfactory solution describing a system of two 
unequal counter-rotating black holes separated by a massless strut, must be depicted by five physical parameters, 
i.e., the Komar masses Mi, M2 of each constituent, their respective Komar angular momenta per unit mass oi, 
a2 (ai and 02 having opposite sign) and the relative coordinate distance R between the centers of the black hole 
horizons. The main difficulty to accomplish this endeavor is the problem of avoiding the NUT sources in order to be 
able to provide the explicit form of cti and (72 in terms of Komar's physical parameters. 

In this work, we first derive a metric describing a system of two rotating black holes by means of the SibgatuUin's 
method [l^ [l^ . Then, we apply the axis condition and solve the corresponding equations for the particular case of 
two unequal counter-rotating black holes separated by a massless strut. We write the solution in terms of ai and 
tT2 , as a 4-parameter subclass of the double-Kerr-NUT problem [ij . Later on, we calculate txi and 02 by using the 
Komar integrals for the masses Mi and the angular momenta Ji. We show that the interaction force between the 
black holes, provided by the strut, has the same form as the one for the static double-Schwarzschild case. Moreover, 
the 5 parameters satisfy an algebraic relation, which generalizes the two statements made by Bonnor [^, [l^ in order 
to remove additional contributions made by the massless spinning rods outside of the sources, i.e., Eq.([T]). In this 
description, the total angular momentum of the system is exactly J — A/iOi + A/2a2. Notice that it contains only the 
contributions from the two sources. 

The outline of the paper is as follows: In Sec. II all the necessary elements to construct a solution for a binary 
systems of Kerr sources are presented. In Sec. Ill a four parametric solution for two counter-rotating black holes 
separated by a massless strut is analyzed. In Sec. IV the parametrization of the solution in terms of the physical 
Komar parameters is accomplished. In Sec. V the concluding remarks are given. 



II. A SOLUTION FOR A BINARY SYSTEM OF KERR SOURCES 



The Papapetrou line element describing stationary axisymmetric spacetimes reads [l9j 

ds^ = /-I [e^'^idp^ + dz^) + p2 V] - f{dt - Lodipf, (2) 

where /, u, 7 are unknown functions depending only on the cylindrical coordinates {p,z). According to the Ernst 
formalism [201, vacuum Einstein field equations for these particular stationary axisymmetric spacetimes read: 

{Re£)A£ = V£ ■ V£, (3) 

where Vand A denote the gradient and Laplace operators, respectively, defined in cylindrical coordinates and acting 
over the complex Ernst potential £ = f + W. For any solution of the equation the metric functions ui and 7 of 
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the line element ([2]) can be obtained from the following system of differential equations 

= -pr'^'fz, uj,=pf-Hp, (4) 

ijp ^ pf-mSp\'-\S.\^), 27,-p/-2Re(fpf,), (5) 

where the bar over a symbol represents the complex conjugate, \x\'^ — xx, and the subscript p or z denotes partial 
differentiation. In order to solve the non-linear Eq. ([3]) , we will use the powerful mathematical technique based on the 
soliton theory known as Sibgatullin's method p^]. The extended double-Kerr-NUT problem jlj can be constructed 
easily by applying this method as it is done in reference (l7j for the case of two bodies, but vanishing electromagnetic 
field ($ = 0). Let us start by defining the Ernst potential on the symmetry axis as follows: 

£ip = 0, z) ^ e(z) = 1 + ^i- + (6) 

The set of complex constants parameters {ej,/3j} are a total of eight real parameters related with the multipolar 
terms. Once we know the value of the Ernst potential on the symmetry axis, the complex potential in the whole space 
can be obtained from the Sibgatullin's integral 

whose unknown function /x(C) satisfies an integral equation 



and a normalization condition 

'-fj^^l, (9) 

» J-i vw^ 



where e{r]) = e{fj) and ^ is representing a principal value integral. In addition, e(^) is the local holomorphic 
continuation of e(z) on the complex plane z + ip, with — z + ip^, rj = z + ipr, V^, t G [—1,1]. Since e{z) is a rational 
function, the corresponding fi{(^) can be assumed to be of polynomial form 
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MC) = Ao + ^A„(e-a„)-i, (10) 

71=1 

where the coefiicients Af) and A„ are determined by the Eqs.®-®, and the constants a„ represent the location of 
the sources on the symmetry axis, they are the roots of the following characteristic equation, see Figll] 

e{z) + e{z)^0. (11) 

Replacing Eq.([6]) into Eq.(fTT|). it is possible to show that the old parameters {ej,/3j} and the new ones {a„,/3j} are 
related through the following relations 

(/3i-/?2)(/3i-^i)(/3i-^2) ' (/32-/?i)(/32-^i)(/32-/32) ■ ^ ' 

After tedious calculations the solution describing the extended double-Kerr-NUT problem can be straightforward 
obtained. The Ernst potential £ and the corresponding metric functions /, uj and 7 can be written in the following 
explicit form [2lj 

E+ _ E+E^ + E+E^ 4Im[^„G] _ E+E^ + E+E_ 

E.^ ^ 2\E.\^ ' " E^E.+E^E.' ' 2\K,mUrJ ^''^ 
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FIG. 1: Location of two unequal Kerr black holes on the symmetry axis represented by the rods of length 2ai and 2(72, where 
the roots satisfy '^an = Q and both bodies are disconnected if the axis condition and R > (t\ + 02 are satisfied. 



These last expressions constitute a solution defined by eight parameters represented by the values of the a„, 
n = 1,2,3,4 and j = 1,2. However, when a„ are set to be real parameters, the solution describes a solution 
for a binary system constituted by two Kerr black holes, where the two horizons are defined on the symmetry 
axis by the intervals ai > z > a2 and > z > a/^. It is important to notice that the above solution was 
constructed assuming asymptotic flatness at spatial infinity, where / — >■ 1, 7 — >■ and w — > (in the absence of 
NUT sources), the metric functions 7 and uj automatically fulfill the following conditions on the symmetry axis: 
7(ai < z < 00) = 7(— 00 < z < Ui) — and uj{ai < z < od) = 0, thus establishing an elementary flatness on the 
upper part of the symmetry axis. 



III. A FOUR-PARAMETRIC SOLUTION 



For the case in which the binary system is located on the symmetry axis in such way that the roots a„ satisfy the 
condition ^ a„ = 0, only seven parameters are needed in order to characterize such solution. In order to get rid of the 
NUT sources between the objects in the lower part of the symmetry axis, i.e., regions < z < a2 and —00 < z < 014, 
thus regularizing the symmetry axis outside the sources, and determining the solution for two counter-rotating black 
holes with a massless strut inbetween, a well-known conical line singularity [ll[, we must impose the following two 
conditions on the metric function uj 

u!{p ~ 0,a2 < z < 03) — 0, u!{p = 0, —00 < z < a^) — 0. (15) 

We note that the second condition in Eqs.([T5|) implies the vanishing of the gravitomagnetic monopole (NUT pa- 
rameter 0), which also can be determined asymptotically by the Ernst potential on the symmetry axis Eq.® as 
follows 



Im[ei + 62] = 0, 



(16) 
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where ei and e-i are defined in Ea. (|T2|) . A straightforward simpUfication performed with the two conditions in Eq.([T5]) 
lead us to the foUowing compact set of algebraic equations 



(17) 



These last two equations reduce the seven parametric solution to a five parametric one and the complete metric can 
be written in a similar form as the one for the double Reissner-Nordstrom problem [2^ [23j . We will restrict our 
solution to a four parametric subclass. Since the solution Eas. (jl3l) - (ll4p involves real constants a„ which determine 
the location of the two Kerr black hole sources on the symmetry axis, we re-parameterize them as follows 
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(18) 



where, as mentioned above, cti and a2 describe the half length of the two rods representing the black hole horizons 
and R is the relative distance between the two centers as it is shown in Figure 1. The lengths ai and CT2 can be written 
in terms of the Komar physical parameters, i.e., the individual masses Mi and M2, the angular momenta per unit 
mass ai and 02, and the coordinate distance R. In our case, the parameters should satisfy an additional relationship 
imposed by Eq. (|17p . The Komar integrals for the individual masses Mi and angular momenta Ji can be calculated 
through the Tomimatsu's formulae [24]: 



M,, 



1 

8^ 



Hi 



uj^^dipdz, = --!-/" w(l + iw*, 



difdz. 



(19) 



The integrals are over the black hole horizons Hi — {a2i < z < a2i-i, < (p < 27r, p 0}, i — 1,2. Moreover, the 
total mass M can be considered as the sum of the individual masses Mi and M2 once the conditions established in 
Eg. psp . for regularizing the symmetry axis outside the sources, are fulfilled. In order to solve the system of Eqs. (ll7l) . 
we make a first order expansion in z~^ of the Ernst potential on the symmetry axis, Eq.(|6]), to calculate the total 
mass M asymptotically, the obtained result reads: 



Re[ei + 62] = -2M. 
Replacing Eg. dT^ into Eg. ([20)1 yields the equation 

Pi+p2 + k+k = -2M, 



(20) 



(21) 



implying several possibilities on the relations between the beta-parameters and the total mass M. The simplest 
choice describing the unequal counter-rotating case is the relation /3i + /32 = —M. A simple calculation lead us to 
the following result 
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£1,2 
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where the subindexes 1 and 2 are associated with + and — signs, respectively. Therefore, writing the Ernst potential 
and metric functions in terms of the parameters M , R, ai and a2 leads to the following four parametric solution for 
two unequal counter-rotating black holes: 
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where 

A = 4CT1CT2 (M^ - ej) (nra + rsr^) + [Af^i?' + e^(i?2 - 2M^)] (n - ra) (rg - r4) + (M^ - 2M^R^ + el) 

X [ei(ri - r2)(r3 - - 2aiCT2(ri + r2)(r3 + Vi)] - 2i5[ai (ri + r2) (ra - r4) - (72 (ri - r2) (ra + r^)], 

r = c7i(Af + e2/Af)[2a2M2(e2-i?')(r3+r4)-i?(M2-e2)'(r3-r4)] 

- (72 (Af - e2/A//) [2aiAf^ (e2 + i?^) (ri + r2) - i?(Af2 + e2f{ri - r2)] 
+ [cri (A/ + e2/Af) (rg - ri) - (J2{M - e2/Af)(ri - r2)] , 

G = 2zr + 4cricr2i?(Af'^-e2)(^i'-2-r3r4)+ai(i?2 + e2)(Af^-e2)^(ri+r2)(r3-r4) (24) 

+ CT2 {R^ - £2) (Ai"2 + e2)^(ri - r2)(r3 + r4) - 2ie25{ri ~ r2)(r3 - r^) - ax (Af + ea/Af) 

X {2(72^? [ea + A//2(Af2 - i?2 _ g^)] (^3 + r^) + [2Af2e^(i?2 _ _^ (-2Af2 - i?2>)g2 _ m^R^\ (r^ - r^) 

+ 18 [R (r3 - r4) - 2^2 (r^ + n)]} + ^2 (A/ - e2/Af) {2a,R [el + M^iM^ - R^ + e^)] (n + '-2) 

- [2Af2ei(i?2 + + (2A/2 _ i?2)^2 _ ^,^4^2j _ _ _ ^ 2(71 in + r2)]}, 

5 := y^(i?2 - M2)(A/2i?2 _ e2)(^^4 _ 2eiAf2 + gZ)^ 
where r„ can be written in the following parameterized form 




ri,2 = 1 U--i?^(7i , r3,i^\ P^+[z + -Rtct2] , (25) 




the indices 1,3 and 2,4 run over + and — signs, respectively. Obviously the solution has not the equatorial 

antisymmetry property in the sense of [isj . the antisymmetry appears only for the case where both constituents are 
equal. 

It is interesting to note that under the transformation 1 <;-)■ 2, z — )■ — z, which exchange the physical properties and 
the position of the constituents, will only change the global sign of the metric function to. The corresponding Ernst 
potential on the symmetry axis now reads: 



e{z) = -, 



-z2 Mz 2A//3 - A/i?2 - 2A/ei T 2e2i? i\/{R^ - M^){M^R^ - el){M^ - 2eiM^ + el) 
e± - 2: T z+ — 2M{MR T £2) ' 



(26) 



The total angular momentum of the system can be calculated asymptotically by means of a second order expansion 
in z~^, it is given by 




_ ^ ; (j?2-M2)(Af4-26iAf2+,2) 
-J - nnr\l APR^ - el ' ^ 

Under the transformation 1 o 2, the total angular momentum changes its sign, i.e., J = — J(i<_>.2). This fact means 
that Eq. (P5|) is indeed a solution for the case of two unequal counter-rotating black holes. On the other hand, from 
the energy-momentum tensor associated with the strut, one obtains the following expression for the interaction force 
between the black holes [HI 

T-'-(e-^"-l)-^^-^^d^ (28) 

where 70 is the constant value of the metric function 7 evaluated on the corresponding region of the strut. It is 
worthwhile to mention, that the interaction force between two identical counter-rotating black holes {Mi — M2 = m, 
ai = —02 = a) in the non-extreme case: M = 2m, (7i = (72 = cr, and in the extreme case: M = 2m, cri = (72 = 0, is 
of the same form [H [3 [11], i.e., 

i?2-4m2- ^^^^ 

Moreover, in the absence of rotation: ai = 02 = 0, ai ~ A/i, (72 = Af2, and M = M-\ + Af2, we recover the 
well-known expression for the interaction force between two Schwarzschild black holes p5l.[27j. i.e., 

AfiAf2 
i?2-(Afi + A/2)2- 
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IV. THE PHYSICAL PARAMETRIZATION 



The relation between the quantities cri , 02 and the physical Komar parameters of the system can be obtained by 
means of the Tomimatsu's formulae Eqs. ((T^ . Let us use the following simplified form of them [2^ [28j: 



M, = ^[^-Ip^C.^a,. -*|p=o,.=c.,._J, J, = y(Mz-aO, « = 1,2, (31) 

for the individual Komar masses and angular momenta, loi are the constant values of the corresponding metric 
function w evaluated over the horizon of each constituent black hole. 

The horizons are defined as null hypersurfaces, i.e., p = 0, ~a\ < z ~ R/2 < ai and p — 0, —(J2 < z + R/2 < a2, 
with a massless strut inbetween. A straightforward calculation leads us to the following system of equations for the 
individual masses and angular momenta of the black holes: 

_ M^ + e2 _ AP-e2 

- -^A^' - ^M-' ^^^^ 



Ji 



Ml / {R + M){MR - £2) (M^ - 2eiM2 + 



M2 {R + M){MR + £2) (M4 - 2eiAP + el) 



2M 



[R^ M){MR + e2) 



2M 



{R- M){MR- €2) 



(33) 



From Eg. (132)) it is easy to see that the total mass M = Mi + M2- However, besides this relation one obtains the 
additional relation 



ol - gI = Ml - Ml 
replacing Ea. (p4)) into Ea. (l33l) leads to the following expressions for ai 
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Eg. (1341) implies the following relation between the five physical parameters: 
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(36) 



This last relation generalizes the two assumptions Eqs.(IT]) made by Bonnor [8l.lT8|. in order to remove the contribution 
arising from the massless spinning rods outside the sources. Notice that Eg. ip?)) accounts exactly for the total angular 
momentum J as the sum of the individual angular momenta of both constituent black holes, i.e.. 



J — A/ifli + M2a2- 



(37) 



As mentioned above, the condition Eq.([36]) means that the contributions to the interaction force made by the two 
semi-infinite massless spinning rods located on the upper and lower part of the symmetry axis have being removed, 
as well as the finite massless spinning rod between the two constituents @, [3]. Hence, Eg . ([281 reduces to the simple 
formula arising for the interaction force between two Schwarzschild black holes. 

Moreover, if the condition Eg. p6p between the five independent parameters is not fulfilled, a proper contribution 
of the spin-spin interaction appears in the expression for the interaction force [28| . 

Notice that we can recover from Eg. ([35]) the case of one isolate black hole by imposing the limit R 00 ox just 
setting to zero the physical properties of the other body, i.e., in this case the corresponding mass. 



TABLE L Particular numerical values for the 4-parameter subclass of the Double-Kerr problem. 
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Table U shows in the first three rows different sets of numerical values for the masses and for the angular momenta 
per unit mass of the black holes. The angular momentum of each component having opposite sign. The fourth row 
displays the case of two equal counter-rotating black holes. The fifth row corresponds to the static case, in which 
the total angular momentum of the system vanishes, i.e., J — 0. In this case, the horizons of the two black holes can 
reach each other and the system evolves into one Schwarzschild black hole. 

Thus, the expressions for £, /, uj and 7, describing our four parametric solution for two unequal counter-rotating 
black holes in terms of physical Komar parameters read: 



where 
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2n21 



= (l/\/2) {R - Ml - M2) Jaj{R + Ml - Ma)^ + al{R - Mi + M2 



(39) 



The Ernst potential on the symmetry axis reads: 
e{z) = — , 



F := 



z^^Mz- 



R-M 
R + M 
R + Ml - M2 



R + M 



R' - M' 



F 



±1 



F^^^^ + ^^alF + alF-^ 



i? - Ml + M2 ' 

(Ti and r„ are given by Eg. ([55]) and Eg.ipSj). respectively. Moreover, oi and 02 fulfill Eg. ([55]) . 



(40) 



A. Thermodynamical properties 

For each component of the binary system, the Smarr formula for the mass [l^l holds, i.e., 

kS 

M, = + 2nia,M, = cr,; + 2nia,M„ i = 1, 2, (41) 
47r 

where Ki is the surface gravity. Si is the area of the horizon, 51^ the angular velocity and the angular momentum 
per unit mass for each constituent black hole. Notice that this last formula implies that Mi > at. In order to calculate 
the values of and Vli , one can use the following relations [28l . i29j 



a,=uj;\ (42) 



being u)i and 7^ the constant values of the corresponding metric functions w and 7 evaluated over the horizon of each 
constituent, while e^'' is negative at the horizon |28l|. By means of the solution Eos. ([55)) - ([M)) . it is straightforward to 
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obtain the following expressions for the angular velocities fli, the surface gravities k^, and the area of the horizons Si 

ai[{R - K'hf - Ml] 



Kl = 



Si 



2A/i(Afi + ai)[{R + M2Y - Ml] ' 
(Ji{R + Ml- M2) 

2Mi{Mi +ai){R + M) ' (43) 
87rA/i(Afi +cri)(i? + Af) 



(i? + Ml ~ M2) 

VL2 — fil(l<^2)7 1^2 — '«1(1^2), 'S'2 — S'i(i<^2)- 

In the limit when the sources are far away from each other, the angular velocities reduce to: 

aiFi 



2M,{M, + a,y 

4 A/2 8Af| 4Af2(A/i2 +3Af|) /I 



4Mi 8M2 4Mi(A/|+3A//f) (I 
P2 ^1-^ + ^2 J^3 + 0[j^, 

notice that the proper contribution to the angular velocity fi^ coming from the angular momentum Ji — MiUi begins 
at the third order of the expansion, i.e., ili ~ 0{\/ R^)) jT2!|. 

Additionally, in the limit Mi = M2 = m, ai = (72 = cr and ai = —02 = a, our solution reduces to the one for the 
case of two identical counter-rotating black holes. The unique a reads 



Therefore, this particular case belongs to a 3-parameter subclass of the double-Kerr solution [l|, where the total 
angular momentum of the system vanishes, i.e., J — Q. 



B. Singularities off the axis 

Since Mi > in ([38|) there exists no ring singularity off the axis. Nevertheless, if one of the masses is negative, even 
when the total ADM mass is positive [30|, the solution Eqs. (l55|) - (I5^ presents such singularities. By setting / = 0, 
this fact can be observed in the following stationary limit surfaces of the Figure [21 




0.00.10.20.30.40.5 0.00.20.40.60.81.0 

P P 



FIG. 2: (a) For positive masses, there exist no singularity off the axis, for the values: ai = 4.973, (J2 = 1.931, M\ — 5, M2 = 2, 
a\ = 1.364, a2 = —3 and ii = 8; (b) If one of the masses is negative, there appears the ring singularity off the axis, for the 
values: a\ = 4.973, 02 = 1.931, M\ = 5, M-i — —2, ai — 0.2, 02 — —3 and R = 8. The ring singularity is located at p ~ 0.69, 
z ~ -4.01. 



The localization of such ring singularity can be calculated as one root of the denominator of the Ernst potential 
Eq.dSHl). 
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V. CONCLUDING REMARKS 



In this work, we present an exact solution describing a binary system constituted by two unequal counter-rotating 
black holes with a massless strut inbetween. We derive a 4-paraineter subclass involving a simple algebraic relation 
between the five physical parameters. This relation generalizes, for systems of unequal black holes, the two assumptions 
made by Bonnor [l^ in order to avoid the contribution from the massless spinning rods outside the black holes. 
Therefore, the interaction force provided by the strut results to be of the Schwarzschild type. This solution reduces 
to the one for the case when both constituents are identical [T3|. 

On the other hand, in the extreme limit: ai ~ and CT2 = 0, the unequal and opposite angular momenta per unit 
mass, in absolute value, are greater than their corresponding positive masses: \ai\ > Mi > 0, i.e.. 



This fact was first pointed out by Herdeiro et al. [3l|, and it can be obtained also as a consequence of the work of 
Varzugin [12]. Additionally, in this particular case, the condition established between the five parameters is satisfied 
only if both constituent black holes are equal, i.e.. Mi = M2 = m and ai = —02 = a. The total angular momentum 
of the system vanishes, and the distance R in which the extremality condition occurs [l4j reads: 

R= \^ \a\>m>0. (47) 

It is worthwhile to mention that expressions Egs. pSl) . relating the masses and the angular momenta, are of the 
same form as the relations presented in [32l| , in the context of a binary system constituted by two extreme Reissner- 
Nordstrom black holes with a strut inbetween. 

The technical details for removing the NUT sources outside the two rotating black holes is not a trivial problem and 
it restricts the possibilities for finding exact solutions to more general problems related to the counter/co-rotating 
cases. 
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